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We construct and investigate a compactified version of the four-dimensional Reissner- 
Nordstrom-NUT solution, generalizing the compactified Schwarzschild black hole that has 
been previously studied by several workers. Our approach to compactification is based 
on dimensional reduction with respect to the stationary Killing vector, resulting in three- 
dimensional gravity coupled to a nonlinear sigma model. Using that the original non- 
compactified solution corresponds to a target space geodesic, the problem can be linearized 
much in the same way as in the case of no electric nor NUT charge. An interesting feature of 
the solution family is that for nonzero electric charge but vanishing NUT charge, the solution 
has a curvature singularity on a torus that surrounds the event horizon, but this singularity 
" is removed when the NUT charge is switched on. We also treat the Schwarzschild case in 

• a more complete way than has been done previously. In particular, the asymptotic solution 

(the Levi-Civita solution with the height coordinate made periodic) has to our knowledge 
£SJ ', only been calculated up to a determination of the mass parameter. The periodic Levi-Civita 

solution contains three essential parameters, however, and the remaining two are explicitly 
calculated here. 
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I. INTRODUCTION 
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String theory predicts that our world has more dimensions than the 4 we see around us every day. The 
usual way out of this apparent discrepancy is to say that the additional dimensions are compactified and 
small. In the simplest case the extra dimensions are just circles. String theory is also a theory of gravity. 
It is therefore a natural and interesting question to ask what is the behavior of the classical solutions of 
General Relativity in spacetimes where one or more dimensions are compactified. This question is however 
surprisingly hard to answer. The periodic analogue of a Schwarzschild black hole in a space of topology 
R3 x SI has been st udied more or less independently by Majumdar pj, Papapetrou 0, Myers 0, Korotkin 
and Nicolai H. Q | as well as by Frolov and Frolov In higher dimensions the full solution is known only 
implicitly [lfllll IH E|. 

The problem becomes even more interesting when one realizes, following Gregory and Laflamme 
that there is a competition between different gravitational configurations with the same mass and charges 
but with different symmetries and even of different horizon topology. Which configuration is stable depends 
on the particular value of the mass and the charges. For instance, a long and thin black string becomes 
unstable and will "decay" into a new configuration with higher entropy. This new configuration could for 
instance be a nontranslationally invariant black string or an infinite array of black holes. The phase diagram 
of black objects in higher dimensional circle space times has a very rich structure [Til ITil \l$L . Only parts 
of it is accessible to analytical methods and the most interesting pices has been calculated only numerically 
[PH liH l22l m 0, or using perturbation theory 0] . It would therefore clearly be of great interest to 
have a better analytical understanding of the various solutions that could appear. 

In this paper we perform a more modest task. We generalize the solution given in to arbitrary 

electric and NUT charge. In doing so we find a new parameter of the solution which has not been recognized 
before. We analyze the physical meaning of this parameter. We also notice that adding NUT charge to the 
solution makes it more well behaved. Maybe this is a feature which would persist in higher dimensions? Our 
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method heavily relies on Weyl coordinates. However, as shown in 26], because of symmetry reasons, Weyl 
coordinates are expected to be useful at most in 4 and 5 dimensions. This is something which would have 
to be overcome if one would like to generalize our results to higher dimensions. 

This paper is organized as follows. Section 2 is a general discussion of properties of the four dimensional 
Einstein-Maxwells equations important for our problem. In Section 3 we discuss the (nonperiodic) Reissner- 
Nordstrom-NUT solution. In Section 4 we find its periodic generalization and discuss the metric in various 
limits. In Section 5 we give our conclusions. 

II. EINSTEIN-MAXWELLS EQUATIONS FOR STATIONARY SPACETIMES 

We write a general stationary metric with timelike Killing vector £ a as 

fjab = -/ 2 MaA*6 + f~ 2 hab, (1) 

where 



/ = V^fa (2) 
Ma = -r^a, (3) 

which implies that h a b£, b = and that /, n a and h a b are stationary fields, i.e. Lie dragged by £ a . We now 
decompose the electromagnetic field strength F a b by defining the electric and magnetic fields E a and B a 
with respect to £ a : 



where *F a b is the dual field strength 



It then follows that 



E a = F ab e (4) 
B a = *F ab Z b , (5) 

*F ab = \e ab cd F cd . (6) 

F ab = 2i?[ Q ^h] — e a b cd B c fid (7) 
*F ab = 2B [al i b] + e a b cd E c ^ d . (8) 

Assuming that F ab is stationary as well as source-free so that *F ab is closed, then both E a and B a are 
gradients of stationary scalars. Hence we introduce electric and magnetic potentials v and u according to 

E a = -V a «, B a = -V„U. (9) 

Moreover, with the electromagnetic field being the only matter source, a twist scalar x, likewise stationary, 
can be introduced byQ 

V aX = e abcd ^y d e + 2(uV a v - vV a u). (10) 

The tensor h ab can be viewed as a metric on the three-dimensional manifold S of Killing orbits and we shall 
use that the full 4D Einstein-Maxwell equations are now equivalent to the 3D equations obtained from the 
action 



S = J d 3 xVh ^R^2h ab lAB D a X A D b X B , (11) 



where ^R, h ah and D a are the Ricci scalar, inverse and Levi-Civita connection of h a b- This action describes 
three-dimensional Euclidean gravity coupled to a certain sigma model. More precisely, the effective matter 
part of the action is that of a wave map from X to a four-dimensional target space with metric 

2 A B id£d£ + if) dipd£ + t/> dipd£ — (£ + £)dijjd'tjj 
da- = "fABaX dX = = — = , (12) 

2 [|(£ + £)- 



where £ and ip are the complex Ernst potentials 

£ = f 2 — v 2 — u 2 + ix, tp — v + iu. (13) 

This is an Einstein metric, whose Ricci tensor obeys 

Rab = -6jab (14) 

If we use as target space coordinates X A the real variables /, A, S, x, where A and S are the amplitude and 
phase of if>, i.e. 

ip = Ae lS , (15) 

the target metric takes the neat form 

V-jsUf-iA'-AW+ te+ZW , (16) 

which makes it evident that it has signature (—,—,+,+). Now, varying the action with respect to h ab 
gives 

^R ab = 2 lAB D a X A D b X B , (17) 

while variation with respect to the target space coordinates X A leads to the wave map equation 

h ab (D a D b X c + T c AB D a X A D b X B ) = 0, (18) 

with T c ab being the Christoffel symbols of jab- For our purposes it is important to note that if the X A 
all depend solely on one free function u>, then eqs. I|17|) and (|18|l become 

dX A dX B 

^R ab = 2 1AB — ^—D a oj D b io (19) 

dw dw 

dX c \, fd 2 X c r dX A dX B \,, , . 

-—h ab D a D b u + — =- + T C AB — — h ab D a cj D b uj - 0. 20 

dw \ du> z dw dw J 

Now, eq. ifflty tells us that the curve X A (u) is a geodesic, since, as h ab is supposed to be a positive-definite 
metric, it cannot happen that h ab D a uo D b u vanishes identically unless lo is constant. Using the freedom to 
reparametrize the geodesic, we can assume that it is affinely parametrized and arrive at the equations 

dX A dX B 

e := JAB—, ; — = constant (21) 

du dui 

d 2 X G | tC dX A _dX^_ Q 

dui 2 did duo 

(3) R ab = 2eD a uj D b oj (23) 
h ab D a D b oj = 0. (24) 

In this paper we focus exclusively on this class of solutions, which we shall refer to as geodesic solutions. 



A. Axisymmetry and Weyl coordinates 

Let us now assume that the stationary spacetime is also axisymmetric with an axisymmetry generator r] a 
that commutes with the timelike Killing vector £ a . We can then introduce coordinates (t, p, z, <f>) such that 
£ a = (d/dt) a , jf = (d/d(f)) a while the three- metric h ab and one- form /x a take the forms 

dl 2 = e 2k (dp 2 + dz 2 ) + W 2 d(f> 2 (25) 
H = dt + nd(p. (26) 
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(29) 



The metric functions f,k,W and f2 obviously depend on p and z only. If we assume that lj also depends 
solely on p and z (although it could have a linear dependence on <f> with constant coefficient, a possibility 
which we do not consider here), the ^-component of eq. Il2.'j[| becomes 

W tPP + W, zz = 0, (27) 

making it possible to choose p and z such that W — p, a choice well-known as Weyl (canonical) coordinates. 
The great advantage of being able to use these coordinates is that the unknown metric function k does not 
enter the Laplace equation l|24|l (nor the more general wave map equation QSJl), which will thus be identical 
to the Laplace equation in flat space, expressed in cylindrical coordinates, for an unknown function that is 
independent of <j>; 

u, pp + p + u), zz = 0. (28) 
With the use of this equation, the remaining components of eq. I|23|) become 

k, p = ep [(w p ) 2 - 
k, z = 2epui tP uJ tZ , 

while eq. (|10|) gives 

Ct n = Ap y PUJ z , . 

S2 iZ = -4p x pw, p , 

with p x being a constant equal to the conserved geodesic momentum associated with the cyclic coordinate 
X of the target space metric. Explicitly 

Px - ^ 4/4 dUJ - (31) 

A recipe for finding a solution of the type considered here is thus to choose an appropriate geodesic of the 
target space and an appropriate solution to the flat space Laplace equation l|28|l . The geodesic then gives 
the four target space functions as functions of to, including the metric component / and the potentials v 
and u which completely determine the electromagnetic field via eqs. (JJJ) and ©. What then remains is to 
integrate the two pairs of equations l|29|l and (|30(l to determine the remaining spacetime metric components 
k and fl. Note that eq. I|28|l is the integrability condition for both pairs of equations (when e 7^ 7^ p x ). 



III. THE REISSNER-NORDSTROM-NUT SOLUTION 

Using Schwarzschild type coordinates, the Reissner-Nordstrom-NUT (RNN) solution with mass M, electric 
charge Q and NUT (gravitational monopole) charge I is 

ds 2 = -f 2 (dt + nd(j>) 2 + f- 2 dr 2 + {r 2 + l 2 )(d9 2 + sm 2 0d(f> 2 ), (32) 

r 2 — I 2 2Qlr 
F = Q-H? j ¥ ^(dt + Q.d(f>) Adr+ ^ sing dd A d<f>. (33) 



where 



(r 2 + / 2 r r 2 + l 2 



a r 2 - 2Mr + Q 2 - I 2 (r-r+)(r-r_) 
J r 2 +l 2 r 2 + p W > 

n = 2lcos9 + Q'. (35) 



Here we have defined 



r± = M ± A, A = y/M 2 + l 2 -Q 2 , (36) 



while fl' is a constant which should be set to —21 (21) to make the half-axis 9 — (9 — 7r) explicitly regular, 
leaving the other half-axis - the Misner string - singular, as d(f> is not a well-behaved one-form at 9 = 0, ir. 
However, it is well-known that both half-axes can be made regular, since changing fl' from —21 to 21 can be 
mimicked by changing the time coordinate from t to t' = t — Al<j). The price to pay is closed timelike curves 
since it requires that t and t' should both be periodic with period 8irl. 
The three-metric h a b for this solution is seen to be 



dl 2 = dr 2 + (r 2 - 2Mr + Q 2 - l 2 )dVt 2 = dr 2 + (r - r+){r - r_)dft 2 



and the three remaining target space scalars are given by 

2l(r — M) 
r 2 + I 2 

Q 



x 

A = 
S 



Vr 2 + 1 2 

arctan | - 
r 



(37) 

(38) 
(39) 
(40) 



Since all the X A depend on r only, this is a geodesic solution. In the non-extremal case A > 0, the geodesic 
is spacelike (e > 0) and explicitly given by 



/ = 

X = 
A = 

S = 



A 



(M sinh lo — A cosh to) 2 + I 2 sinh 2 cj 

I A sinh 2u> 
(M sinh Co — A cosh a)) 2 + I 2 sinh 2 a) 
Q sinh<I> 

(Msinhw — A cosh Co) 2 + i 2 sinh 2 tj 
I sinhw 



2lA(e 2Q 



l 2 (e 



>(e a 



arctan 



V M sinh lo — A cosh co 



arctan 



l 2 (& 



r + e" 



(41) 

(42) 
(43) 

(44) 



where oj is the arclength parameter for the geodesic, and here the function that satisfies the Laplace equation 
(PHl . namely 



arccoth ( 



M 



\ A 



However, we shall instead think of the geodesic as parametrized by the rescaled affine parameter 

M _ 

LO = — LO, 



(45) 



(46) 



which implies that the norm of the geodesic tangent vector has norm e = A/M. The reason for this is that 
we will then straightforwardly be able to treat the extremal case - which corresponds to a null geodesic 
(e = 0) — as the A — > limit of the non-extremal case. Indeed, taking this limit for the above formulae, we 
find that the spacelike geodesic goes over into a lightlike one; 

/ 



X z 
A 



^M 2 {1 


- LO) 2 + l 2 L0 2 




2lMuo 


M 2 (l 


- LO) 2 + l 2 L0 2 




Quo 



S — arctan 



^M 2 (1-lo) 2 + Plo 2 
Ilo 



M(l - lo) 



(47) 
(48) 
(49) 

(50) 



G 



while u goes over into 

M 



-M 



(51) 



This limiting procedure would not have worked, had we used uj as the affine parameter. Moreover, u has 
the large r asymptotic behaviour 

cj = -— + 0{r- 2 ), (52) 
r 

which means that it is u>, rather than Q, that corresponds to a Newtonian gravitational potential. However, 
although we will think of the geodesic as parametrized by uj, we will often work with Q in subsequent 
calculations, since the factor M/A would otherwise often appear merely as an annoying appendage. 

For future reference, we here finally calculate the conserved geodesic momentum p x according to the 
formula Ij31|l . to find the simple relation 

Px = ~. (53) 



A. Charges 

A special and somewhat surprising feature of the RNN solution is that although no sources are included 
for the electromagnetic field, the result of calculating the electric charge Q c \ by integrating *F a b over a two- 
sphere of constant t and r and dividing by A-k is not what one would naively expect from Stoke's theorem, 
i.e. the result is not independent of the choice of two-sphere but depends on the radius; 

r 2 -I 2 

Qel = Q^-~j2- (54) 

We also calculate the magnetic charge <5 mag by integrating F a i, over the same sphere to find 

2Qlr 



! mag 



r 2 + l 2 



(55) 



Obviously we have Q e \ = Q, Q mag = asymptotically as r — > oo, so at infinity the solution is purely electric 
with charge Q. It is worth noting that the relation 



?! + Qmag = Q 2 (56) 



holds as an identity, for all r. The reason why it is possible that electric and magnetic charges depend on 
r is that the gravitomagnetic vector potential £ld(j) enters not only the metric but also the electromagnetic 
field strength (J2HJl. Thus the latter is singular on the Misner string, i.e. on 8 = or 9 = ir (or both) 
depending on how Q' is chosen. As the Misner string goes through every two-sphere that we integrate over, 
we cannot expect Stoke's theorem to imply that the charges should be independent of the choice of sphere. 
Interestingly, if we refrain from removing the singularity by introducing periodic time, the Misner string 
must be thought of as a concrete physical object as it carries electric and magnetic charge density. 

B. Transformation to Weyl coordinates 

For the RNN solution, the transformation from Schwarzschild type coordinates to Weyl coordinates reads 



p = x/{r - M) 2 - A 2 sin 9 (57) 
z = (r- Af)cos<9, (58) 
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with inverse transformation 



r = A + M (59) 



z 



cos9=-, (60) 
A 



where 



1 



A=-(A++A_), A± - y/p 2 + (z±A) 2 . (61) 



The functions u>, k and Q are now given by 

M , /A- AN 

"=2A 1 *{xTa) (62) 
^iln(^-f), ,= i(A + -A_) (63) 

n = ^ + n'. (64) 

A 

In the non-extremal case A > 0, the function lu is the potential of an infinitely thin rod located at p = 0, 
\z\ < A and having a line density M/(2A) per unit length, which means that the line density can be 
interpreted as the mass line density. Note that in the extremal case A = 0, the function k vanishes identically 
leaving the three-metric h a b hat. Moreover A reduces to y/ p 2 + z 2 and u> to the standard monopole solution 
—M/X. 



IV. PERIODIC ANALOGUE OF THE REISSNER-NORDSTROM-NUT SOLUTION 

We will now construct compactified versions of the RNN solution, which generalize the compactified 
Schwarzschild solution discussed in EJ0,0, E|. The approach we take can be summerized as follows: 

• For every member of the RNN family of solutions, compactify the function to by using canonical Weyl 
coordinates and taking the z-periodic analogue of the original solution to the Laplace equation l|28|) . 
This is straightforward due to the linearity of the latter. 

• Insert u> into the same geodesic that defines the original RNN solution. 

• Integrate eqs. Ij29(l and (|30|l to find the remaining metric functions k and f2. 

Before we start, some general remarks about the first of these steps are in order. As proved by Korotkin 
and Nicolai 0,13, if ° w (p, z ) is a solution to the Laplace equation with the asymptotic behaviour 

%{p,z) = -— + 0{f- 2 ) asr^oo, (65) 
r 



where f = y/ p 2 + z 2 and M is some constant (in our case the mass), then the series 

{M 
—— if n ^ 
L\n\ , (66) 

if n = 

is convergent for any (p, z) such that (p, z + nL) does not coincide with a singular point of °uj{p, z) for any 
integer n. The resulting function uj then obviously defines a z-periodic solution to the Laplace equation (|28J) 
with period L. The constants a n are essential as they make the series converge, but one cannot say that 
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they are uniquely determined as one could add to a n any other n-dependent constant b n that falls off faster 



than 



as n — > oo so that 



B:= £ b n 



(67) 



is a finite constant. In other words, the z-periodic analogue of the function ui(p, z) is a priori only defined 
up to an additive constant B. Korotkin and Nicolai fix this constant by taking eq. (|66[) as it stands (B = 0) 
to define lu(p, z), but this choice has not been made by the other workers that have studied the compactified 
Schwarzschild black hole; while MyersQ makes an explicitly different choice of a n leading to a certain non- 
zero B, Frolov and FrolovQ instead use a Green's function method to compactify which in effect corresponds 
to a third choice of B. Now, since the function w(p, z) is periodic in z, it does not make sense to study its 
behaviour for large f , but for large p it has the asymptotic behaviour 



2M 



lnp + (9(l). 



(68) 



If uj would tend to a constant in this limit, it would be natural to make that constant vanish by a suitable 
choice of B. However, since u> instead diverges logarithmically (unless M — 0), we see no physical motivation 
for fixing B at any particular value and hence we will keep it as a free parameter. 

We will now implement the compactification scheme as outlined above. We thus take °u> to be the function 
given by eq. (|62|l . According to the above, its periodic analogue is 



E 



2A V A„ + A 



B, 



(69) 



where 



+ (z„ + A)2 + + ¥n - A) 2 

M ifn^O 



L\nY 
if n = 0. 



z n = z — nL 



(70) 
(71) 



We take z to have the range z £ [-L/2, L/2] with the end points of the interval identified. The horizon, 
where uj diverges to minus infinity, will thus be located at p = 0, |^| < A, just as in the non-compactified 
case. Since we shall not consider the case when the horizon overlaps itself, we require that A < L/2. It will 
prove useful to introduce the dimensionless variable 



2A 



(72) 



having the range (3 S [0, 1]. For convenience, we will, for the time being, work with uj and use the dimen- 
sionless coordinates 



A' 



y 



z 

A' 



(73) 



Clearly uj is then given by 



E 



1 h/^ll 

2 V A' + 1 



B, 



(74) 
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where 



a: =~ 



M 



^x 2 + (y n + l) 2 + ^x 2 + {y n - 1) 2 J , y n = y-2nf3- 1 
2\n\ 

if n = 0, 



(75) 
(76) 
(77) 



We shall now explicitly evaluate u> close to the symmetry axis x = 0. As the behaviour on the horizon section 
|y| < 1 of the axis is different from the off-horizon section 1 < \y\ < /3 _1 , these two sections will have to be 
treated differently, but in both cases we shall express the result in terms of the function 



having the properties 



<P/3(t)<Pl3(-0 = - 



1 



(78) 

(79) 
(80) 



Moreover, for reasons that will become clear, we shall replace the free constant B by a free constant u 
according to 



B = u-P 1 - l -\n [2^(2)], 
where 7 is the Euler constant. Now, for small x and \y\ < 1, we find that Q is given by 



„2« 



AH{y) 



x 2 + 0(x 4 ), 



where H (y) is the function 



H(y) 



2^(2) 



ipp(l +y)tpp{\ - y)' 
For 1 < \y\ < P -1 , on the other hand, one finds that 



,,2 a 



2^(2) <pp{\y\-\) 



+ 0{x 2 ). 



(81) 



(82) 



(83) 



(84) 



This is all that is needed to generate the whole solution. Indeed, for both axis sectors, cu is a solution to the 
Laplace equation which can be expanded in terms of a constant C and a function A(y) as 



(-l) fe d 2k A(y) x2k 



k=0 



2 k k\) 2 dy 



2/, 



(85) 



Inserting this expansion into the right hand sides of eqs. i|29[l and (|30|) (using eq. (|53|l for the latter), we find 
that k and Q are determined up to additive constants k' and f2'; 



k = k' + C 2 lnx + 2CA{y) - - [CA"{y) + A'(y) 2 ] x 2 + 0(x 4 ) 



n = n' + 21 



2k(-l) k d 2k - 1 A(y) x2k 



^ {2 k kl) 2 dy 



2k-l 



(86) 
(87) 
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where we refrain from giving the whole series for k since it depends nonlinear ly on u>, although it would 
not be difficult to write down at least a few more terms. The value of k' has to be determined from the 
requirement that the axis be regular. Starting with the off-horizon section, regularity means absence of a 
conical singularity which is in turn means that k should vanish as x — ► 0, at least if we assume that Vt 
vanishes in that limit, which can always be achieved as we shall see below. Thus, since the constant C is 
zero in this case, we directly obtain k' = 0. For the horizon section we instead have C = 1. In this case we 
may determine k' by inspecting the horizon metric, which reads 

ds| = R 2 

while R is the constant 

R = sfi% + P e~ a (89) 

Now, since 

JJ'(±1)=T2, (90) 

it follows that for the horizon's polar points y — ±1 to be regular, we must set 

fc' = -u + ln4. (91) 

When we are on the horizon \y\ < 1, the constant f2' in eq. I(87|l precisely corresponds to the constant fi' 
in eqs. I|35|l and (|64f> and should be set to —21 (21) to make fl vanish at the horizon pole y = 1 (y = — 1). 
Continuity at y = ±1 then requires us to set Q' = for y > 1 (y < 1) and f2' = —41 (Of = 41) for y < 1 
(y > 1). Again, these two choices do not correspond to different physics, but are related by a change of 
time t —* t' = t — 4l(j) which results in time periodicity with the period 8ttL However, note that unlike 
the non-periodic case, both choices imply that fl has a jump discontinuity where y = fi~ x is periodically 
identified with y — —(i~ x . This discontinuity is of course not of a physical nature either, since we can shift 
it to any other position on the y-axis. For instance, we may put the discontinuity somewhere on the horizon 
and obtain a completely regular off-horizon section on which f2 is everywhere zero. 

A. Fourier expansions 

In principle, the full RNN solution is given as the x = expansions H85JI - i|87|l . but since these series 
converge very slowly for large x and are quite useless when it comes to determining the asymptotic behaviour 
as x — > oo, we shall now follow Frolov and Frolov and represent the solution in terms of fourier series, thus 
making explicit use of the periodicity of the coordiante y. Still using the rescaled coordinates (x,y), we 
assume fourier expansions for the functions u> and k of the forms 

oo 

u> — uj q (x) cos(irq(3y) (92) 

q=0 

oo 

k = k q (x) cos(-Kqj3y) (93) 

Since cu is a solution to the Laplace equation 1)28(1 . it follows that 

uj' q \x) + x- 1 ^) - (irqP) 2 Q q (x) = 0. (94) 

Disqualifying solutions that diverges exponentially as x — > oo, we obtain that H)q(x) is a linear function of 
In x while, for q > 0, 



2(k'+u) j 2 



u) q (x) (x K (nq(3x), 



(95) 



where Kq is the modified Bcsscl function of the second kind. For small x, we have 

2 \ 



-K (nq(3x)= 1 + ln(^^j +0(x 2 



We can now use that we already know u> for small x to find that 

wo (a:) = u- ~ In [2^(2)] +0 lnf-^ 



Wg(x) 



7TQ 



sin(7rg/3) K (nq(3x) for g > 0. 
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(96) 

(97) 
(98) 



We use this fourier expansion for Co to plot the metric function / for two choices of the parameters M , Q 
and Z (figure nj. 

Due to the nonlinearity of eqs. I|29|) which governs fc, we shall not attempt to give its full fourier series. 
However, since all k q (x) with q > decay exponentially as x — > oo, it will suffice to determine fco(x) to 
obtain the asymptotic behaviour in that limit. Now, the zeroth order fourier term of the first of eqs. (|29|l 
reads 



k' Q (x) = x [(^) 2 - {Co, z )\ = x \co' Q (x) 2 + ~ £ [^(x) 2 - (tt?/?) 2 ^^) 2 ] 



9=1 



Inserting the determined expressions for the u) q (x) gives 

f 1 °° r 

fe (x) = ,3 2 x <^ — + 2^ sin 2 (7rg/3) Ki^qfay - K (irqf3xY 

[ x q =i 1 ) 

which can be integrated to yield 

{oo 
In x + x 2 sin 2 (nqp) 2K X (<!rq/3x) 2 - K {nqf3x)K 2 (irq(3x) - K {nq[3xY 
9=1 



(99) 



(100) 



(101) 



where is some constant yet to be determined. Using now that 

2 



2K 1 (nkPx) 2 - K (nkPx)K 2 (TTk(3x) - K (nk(3xf 



(tt/c/3) : 



1 + 7 + In 



irk/3x 



+ 0(x 2 ), (102) 



we can compare ko (x) as given by eq. I|101|) to the zeroth order fourier term of the exact expression for k 
we have near x — to find, after some manipulations, that fc* is the /3-dependent constant 



K = (3 2 



1 + ln 



-/?ln[^(2)] + / ln[^(2)]d/3'. (103) 
Jo 

Finally, from eqs. H30(l and (|53|l we obtain that O is given by the following "quasi-fourier" series: 

oo ^ 

= 0' + 2Z,3y + M V — sin(7rfc/3) xi^i (7rfc/3x) sin(7rfc/3y), (104) 



k=l 



where O' is the same constant that appears in eq. (f%T|) for the case \y\ < 1. 



B. Properties of the event horizon 

As noted above, the (outer) event horizon has the geometry 

dy 2 



dsfi = TV 



H(y) 



H(yW 



TZ = e- u \/r 2 + + P. 



(105) 
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FIG. 1: The metric function / = y/—gu, plotted for M = 0.15L, Q = I = as well as for 
M = Q = I = 0.15L, with u — in both cases. As the two top plots show, the behaviour 
of / for the two cases is very similar close to the axis p = 0. On the other hand / behaves 
qualitatively different when a larger range of p is included, as indicated by the plots below; 
in the first case / increases monotonically with p towards infinity, while in the second case / 
reaches a maximum and subsequently decreases towards zero. The latter behaviour is generic 
as long as I 7^ 0, regardless of Q, but with I = 7^ Q, the maximum is replaced by a divergence, 
as discussed in subsection II V Dl 



It directly follows that the horizon area and Gaussian curvature (half the Ricci scalar) are given by the 
simple formulae 



A = 4+K 2 

K = - l -H"{y)U- 2 . 



It is worth noting that H(y) behaves for small as 



H(y) = l-y 2 - -C(m ~ V 2 ) 2 P 3 + 0(/3 5 ), 



while for (3 close to unity it behaves as 
H(y)= 4(1-/3) +4 



*ii±uH*ri^ i+2 7+ i 



(1 - /3) 2 + 0((1 - P) 3 ), 



(106) 
(107) 

(108) 

(109) 



where $ is the digamma function (the logarithmic derivative of the gamma function). Clearly, for p = 0, the 
metric is an exact two-sphere of radius 1Z. Note that (3 — ► corresponds to two independent and physically 
different limits, namely the extremal limit A — > as well as the limit of infinite coordinate diameter L — > 00. 
This is natural since both limits, makes the proper distance between the poles of the black hole infinite, 
keeping the black hole from being distorted "by itself", i.e. by the gravitational field from its periodic 
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copies. As (3 is increased, the sphere becomes deformed in a prolate manner. For small (3 the deformation 
is very small, however, as the first correction term in eq. (|108|1 enters at order /3 3 . As (3 is further increased 
the deformed sphere gets stretched out to a long cigar, and for /3 close to unity the geometry is almost 
everywhere a flat cylinder of length TZ/^l — j3 and circumference Airll^/l — f3, except at the poles y = ±1 
where the curvature is very high (note that the mean Gaussian curvature is always TZ~ 2 ). Up to a constant 
conformal factor the geometry only depends on the parameter /3 and hence we can refer the reader to Frolov 
and Frolov 6J for a further discussion of its properties, as well as embedding diagrams. Now, the horizon area 
A is clearly modified by a factor e~ 2u compared to the non-compactified case. This is in agreement with the 
general treatment of distorted black holes with (electric) charge as presented by Fairhurst and Krishnan 27], 
based on Geroch and Hartle's treatment of the case without charge[28|. Although it may not have been 
obvious from the way we introduced it, the parameter u (denoted the same way in [2jj but with an unbarred 
u in and Q) can be defined as 

u = Slu\ p=0 ,z=±a, Su) = Q — °Q, (110) 

i.e. u is the value of 8uj evaluated at either of the poles of the event horizon, with 5uj being the difference 
between the function uj for the deformed black hole (whether the deformation is due to compactification as 
in our case, or to external matter) and the function °w for the undeformed RNN solution with the same M, 
Q and /. Since asu-* — oo, we have 

/=— ^^ + 0(e 3 *), (111) 

it follows that u is closely related to the change in horizon pole redshift factor that the deformation produces. 
In the Schwarzschild case Q = I = 0, the relation is more direct since / = e ffl in that case. The term "horizon 
redshift factor" must not be taken too literally, however, since it normally refers to an asymptotically flat 
situation when the value of / _1 gives the redshift factor between the point of evaluation and infinity where 
/ = 1. Anyhow, we shall refer to u as the redshift parameter. In tabled we display how u has in effect been 
fixed by previous workers that have studied the compactified Schwarzschild solution. As we see no reason 
why any one of these choices should be better than the other two, we prefer to keep u as a free parameter. 

Other quantities that are interesting to calculate on the horizon are the surface gravity n and Komar mass 
AfKomar with respect to the Killing vector £ a , as well as the electric and magnetic charges and potentials. 
We find that these quantities are given by 



A 

^TZ 2 


(112) 


-^Komar = A 


(113) 


r\ - I 2 

Qc\-Q r 2 ++l 2 


(114) 


2Qlr+ 
Vmag ~ r\ + P 


(115) 


Qr+ 
V r\ + l 2 


(116) 


Ql 

U ~ r 2 + I 2 


(117) 



We note that the surface gravity is constant over the event horizon and thus the zeroth law of thermodynamics 
holds. Moreover the electric and magnetic potentials are constant as well. In fact the electromagnetic 
charges and potentials take the exact same values as in the noncompactified case. For a discussion of 
further thermodynamic properties, we refer the reader to the general framwork for distorted charged black 
holes I27], Hi| . Unfortunately, to our knowledge no such framework exists in the case of nonzero NUT charge, 
but providing it here would be beyond the scope of this paper. 
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WUI j\.c±. o 




Myers [£] 


±ln[2 W (2)]+ln 


Korotkin and Nicolai [4, 5^ 


|ln[2^(2)]+/3 7 


Frolov and Frolov \6j 


iln[2^ (3 (2)]+/31n(47r) 



TABLE I: The redshift parameter w as used in previous studies of the compactified Schwarzschild solution. 



C. Proper distance between black hole poles 

A quantity that characterizes the off-horizon section of the symmetry axis is the proper spatial separation 
between the poles of the black hole. It can be calculated as 

L scp = 2Aj f-Hy, (118) 

with / evaluated at x = 0. To get a general idea of how this quantity depends on the choice of charges 
M, Q and I, we have held the quotients Q/M and l/M fixed at four different values and, with the redshift 
parameter u set to zero, plotted L sep /L as a function of j3 (figure EJ. Since L scp depends nontrivially on 
how one chooses u to depend on (3, ourplot for the case Q = I = looks different than the one given by 
Frolov and Frolov In particular, in |(| L sop tends to a finite value rather than zero as (3 — > 1, which is 
counter-intuitive since in that limit the polar points reach each other as the event horizon fills the whole 
axis. Again, this is directly related to the behaviour of u in the same limit. However, any choice of u which 
stays finite results in a vanishing L sop as /3 — * 1, just like for u = 0. To show the effect of different choices 
of u, we have also plotted L sop as a function of j3 for the different choices of u collected in table |U with the 
choice u = again included for comparison (figure |3j) . 



D. Large distance asymptotics 

In the limit p — > oo, our compactified black hole solution approaches an exact solution corresponding to the 
field of a homogeneous line mass with electric as well as NUT line charge. This solution is obtained simply 
by truncating all terms that fall off exponentially with p, which is easy to do given the fourier expansions 
derived in subsection IIVAI Explicitly, 

£ D = n-Il 11 [2^(2)]+/31n(^) (119) 

k = K +/3 2 ln(-|) (120) 
Alz 

n = — + n', (121) 

ij 

so the spacetime metric for our compactified Reisser-Nordstrom-NUT solution, in the nonextremal case 
A > 0, rapidly approaches 



ds 2 = -f 2 (dt + nd(/>) 2 + f-' 2 C 2 (p/ Po ) 2p (dp 2 + dz 2 ) + p 2 ^ 2 

where 



(122) 



f = (r± LzL (123) 

[r.(p/p Q )0 - r + {p/p )-P] 2 + P [{p/p o )0 - (p/po)-?} 2 



15 



2.5 



2 - 



Q = 0.75*M, 1 = 



Q = M, l = M 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

P 

FIG. 2: Proper distance between black hole poles as a function of the parameter f3, for u = 
and four particular choices of the quotients Q/M and l/M. 



As before, r± = M ± A with A = \J M 2 + I 2 — Q 2 , but we have here also introduced the constants 

Po = [2^(2) e- 2a p2L 



C = 



(fp(2)e 
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If we set both the electric charge Q and NUT charge I to zero, the function f 2 simplifies to 



(124) 
(125) 

(126) 



which shows that this is a special case of Levi-Civita's static cylindrically symmetric vacuum solution. As 
discussed for instance by Bicak et al. |29j |. the Levi-Civita metric contains two essential parameters m and 
C and can be written in the standard form 



ds 2 = -p 2m di 2 + p 2m ^ m - l \dp 2 + dz 2 ) + i-^ 1 -™)^ 2 

C 2 



(127) 



where the coordinates t, p and i as well as the conicity parameter C in general carry dimension in a 
nonstandard way. The conicity parameter can obviously be set to unity by rescaling the coordinate <j) as 
4> — > C(j), but since we insist that ^ be a periodic coordinate with the standard range 4> € [0, 2</>), the value of 
C cannot be viewed as a coordinate choice. In fact, since we are dealing with a compactified spacetime for 
which z is also a periodic coordinate with a fixed range, the general static cylindrically symmetric solution 
- what one might call the compactified Levi-Civita solution - has instead three essential parameters, which 
for instance can be taken to be our j3, po and C, with j3 directly corresponding to m. 



16 



1.8 



1.6 



1.4 



1.2 



L / L 

sep 



0.8 



0.6 



0.4 



0.2 





I I I I I 

Korotkin and Nicolai 


- 


. ~ — — Myers " 




N^ero redshift 




^parameter 


- 


~n Frolov and \ 




~- Frolov \ 







0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



0.7 



0.8 



0.9 



P 



FIG. 3: Proper distance between black hole poles as a function of the parameter 0, for 
Q = I = and four different choices of the redshift parameter u, three of which being taken 
from the previous studies of the compactified Schwarzschild solution. 



Setting the electric charge Q to zero while having a non-vanishing NUT charge I makes it possible to 
reexpress / 2 as 



2A 



1 



r + (p/p )^-r-(p/po)- 2 P I cosh [2(3 \n(p/ Pl )Y 
where (note that r_ < in this case) 



(128) 



/'I = ( — | Po 



(129) 



As should be expected, this vacuum solution is the "cylindrical analogue of NUT space" that was constructed 
and studied by Nouri-Zonoz 32J. According to Nouri-Zonoz, the "gravitomagnetic charge per unit length" is 
(up to sign) equal to the coefficient of z in the expression for fl, i.e. in our case Al/L. However, what we have 
constructed here is the asymptotic field of a black hole with NUT charge I in a space with a compactified 
dimension of (z-coordinate) length L, so we expect the NUT charge per unit z-length to be l/L. We believe 
that the discrepancy of a factor of 4 is only a matter of definitions (i.e. Nouri-Zonoz' "gravitomagnetic 
charge" is four times NUT charge, when the latter is definied so that its value for the RNN metric Ij32|l is I. 

Setting I = but keeping a non-zero electric charge Q, the general asymptotic solution i|122|) instead 
becomes the cylindrically symmetric Einstein-Maxwell solution with a purely electric Maxwell field. This is 
one of the solutions constructed and studied by Richterek et al. j^lj as well as Miguelote et al. ^31]. The 
function f 2 can in this case be written as 



sinh^ln(p/p 2 )] 



Pi 



r \2f> 
— ) /*'• 



(130) 
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Note that f 2 diverges at a finite cylinder radius p = p?, indicating that there should be a curvature singularity 
there. This is indeed the case, as noticed in [301 l31| . but let us see this explicitly by calculating the 
energy density of the electromagnetic field with respect to the unit timelike vector field u a = where 
£ a = (d/dt) a , as before. The result, in the general case with all three charges M, Q and / arbitrary, is 



Mem = u a u b T ab = \u a u b G ab = 2 



Cp Lj 



(p/Po) 



-2(l+/3 2 ) ji 



(131) 



Clearly /jem diverges at finite p exactly when f 2 does. However, the curvature singularity at p = p2 
disappears as soon as a NUT charge I is turned on, since the denominator of the right hand side of eq. I|123(l 
is then the sum of two squares which cannot vanish simultaneously. Of course, one can make the curvature 
grow arbitrarily large close to p — pi by choosing the quotient l/Q sufficiently small. This feature is not only 
present in the asymptotic solution discussed here, but also in the full compactified RNN solution with the 
z-dependent fourier terms turned back on. The reason for this is that the range of the function ui changes 
from oj £ (— oo,0) to oj e (—00,00) when we compactify the spacetime by making the coordinate z periodic; 
In both cases u> — > —00 occurs at the event horizon, i.e. at p — 0, \z\ < A, but p — > 00 sends u> to +00 rather 
than zero in the compactified case. In particular, for Q ^ = I, the compactified u> will take on the positive 
value 



^crit 



2 \ r - 



(132) 



for which f 2 diverges, as can be seen from eq. ()41(l . Again, this gives a curvature singularity which however 
is regularized by turning on a NUT charge I, however small. The difference, compared to the asymptotic 
solution, is that the curvature singularity now occurs on some torus p = f(z) which is not a simple flat torus 
p = constant. The asymptotic solution also has a curvature singularity on the axis p — for j3 S (0, 1), 
but this is in contrast to the full compactified RNN solution which by construction is regular on that axis. 
Another interesting feature related to the NUT parameter I is that although the spacetime metric 1)122(1 is 
cylindrically symmetric in the sense that no physically measurable quantities depend on the z, this coordiate 
still appears in the expression 1(121(1 for f2. What is interesting is that if we demand that the metric be 
explicitly periodic in z with period L, then the metric must be unchanged under U,' — > tt' + Anl for any 
integer n. This is true, granted that t is periodic with period 87rZ (just like for the non-compactified RNN 
solution we started out with), since we can mimick the shift in Q! by letting t — > t + Anl(p. 

Turning finally to the extremal case A = 0, we must use u> = Mu>/ A instead of Co before taking the limit 
A — > 0. The limit then gives 



2M 
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u' + 7 + lnf — ] , u' — lim 3 1 u 
V2L/J 0^0 



k = 

ilz 

Li 

so the asymptotic spacetime metric is in this case 

ds 2 = -f(dt + ndcf>) 2 + f- 2 (dp 2 + dz 2 + p 2 d<p 2 ) 

with, using eq. l|17| . 

1 



1 - ^- Mp/p3) 



Mp/ps) 



p 3 = 2Le~ u -\ 



(133) 
(134) 
(135) 

(136) 
(137) 



Just like in the nonextremal case, this asymptotic solution fails to be asymptotically flat, albeit in a weaker 
manner. Also, the discussion concerning the curvature singularity at finite p applies here as well. In 
particular, from eq. 147(1 it follows that if I = 0, the asymptotic solution has a curvature singularity at 



p = p3e L/(2M)^ 



but the singularity is removed when the NUT charge I is turned back on. 
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V. 



CONCLUSIONS 



Having constructed and analyzed what we consider to be the natural periodic analogue of the Reissner- 
Nordstrom-NUT solution, a question that naturally arises is whether or not any given asymptotically flat 
black hole has a periodic analogue and, if so, whether or not the periodic analogue is unique, given the period 
L of the compactified coordinate. By identifying the target space geodesic that the non-compactificd solution 
defines and assuming that the very same geodesic should be used for the compactified version, we have found 
that the problem essentially reduces to finding the appropriate axisymmetric solution to the Laplace equation 
in flat space, just as in the Schwarzschild case. There is thus a simple superposition principle at play, which 
makes compactification tractable when Weyl coordinates are used. However, the low dimension requires a 
regularization procedure which introduces a one-parameter ambiguity as there is no obvious way to fix the 
freedom to add an arbitrary constant to the solution to the Laplace equation. In all previous works on the 
periodic Schwarzschild solution, this constant - the redshift parameter - was always fixed at an early stage 
and without any physical motivation for the specific choice made. As we have seen there are three different 
choices in the existing papers on the subject, which we take as support for instead keeping the parameter 
free, or at least fixing it at a later stage from physical considerations. One way of fixing it would be to 
require that a set of relations that hold for horizon quantities (for instance the relation between the horizon 
area and the black hole charges) should remain the same after compactification. If one could motivate why 
this should be the case, our method would give a neat one-to-one correspondance between standard and 
periodic Reissner-Nordstrom-NUT black holes. Of course, we have not given any proof that our method 
is correct, in some suitably defined sense. One could perhaps argue that it is not, based on the fact that 
our solution family has a curvature singularity surrounding the horizon in the case of electric but no NUT 
charge. However, we see no way of avoiding such a singularity. Indeed, the singularity in question is actually 
a general feature not just for our full solution with a nontrivial dependence on the periodic coordinate, 
but also for the cylindrically symmetric solution that is approached for large cylinder radii. Since there 
simply is no other cylindrically symmetric solution available that could work as an asymptotic solution 
for a compactified Reissner-Nordstrom black hole, it seems at least likely that the exponentially decaying 
corrections to the asymptotic solution cannot smooth out the singularity. 

One can think of several ways to extend our work while staying in four dimensions. One way would be to 
include a dilaton field with arbitrary coupling parameter. This would be straightforward to do, as our method 
would still be applicable with the dilaton simply entering as an extra target space coordinate. It would be 
interesting to see whether the inclusion of a dilaton could prevent the above-mentioned singularity from 
occuring. A general framework for distorted charged dilaton black holes has been provided by Yazadiiev[33|. 




It would be even more interesting if the periodic analogue of the Kerr solution could be constructed. Some 
indications of how this could be done were given in 5], but no explicit formulae were given there, except for 
a few first steps. Since the noncompactifed Kerr solution traces out a target space two-surface rather than a 
geodesic, it would no longer be possible to linearize the problem and thus more sophisticated methods, such 
as Backlund transformations, are needed. 

It would of course be even more interesting if our work could give some insights into how things work 
in higher dimensions. Consider the Einstein-Maxwell equations for a static (d + l)-dimensional spacetime. 
Dimcnsionally reducing with respect to the static Killing vector, one arrives at Einstein gravity in d euclidean 
dimensions, coupled to a two-dimensional (three-dimensional if a dilaton is included) sigma model. One 
question one could ask is whether or not the target space geodesic that the (d + l)-dimensional Reissner- 
Nordstrom solution corresponds to would also be traced out when the solution is compactified on a circle. 
The work of Harmark and Obers actually suggests that this could be the case, as their ansatz is such that 
the target space coordinates are functionally dependent, thus corresponding to geodesic solutions. This does 
not mean that it would be easy to make analytical progress for d > 3, but nevertheless it could give a 
coordinate invariant way of thinking about the ansatz of these authors which in turn could lead to a deeper 
understanding of the general problem. It would be very interesting to investigate this matter further. 




19 



Acknowledgements 

The work of RvlJ was supported by the Ministry of Education of the Czech Republic under the project 
MSM 0021622409. The work of MK was supported by a post doctoral fellowship of the Faculty of Science 
of Masaryk University. 



S. D. Majumdar, "A Class Of Exact Solutions Of Einstein's Field Equations," Phys. Rev. 72 (1947) 390. 
A. Papapetrou, Proc. R. Irish Acad. A51, (1947) 191. 

R. C. Myers, "Higher Dimensional Black Holes In Compactified Space- Times," Phys. R ev. D 35 (1987) 455. 
D. Korotkin and H. Nicolai, "A Periodic analog of the Schwarzschild solution," arXiv:gr-qc/9403029 
D. Korotkin and H. Nicolai, "The Ernst equation on a Riemann surface," Nucl. Phys. B 429 (1994) 229 
|arXiv:gr- qc /9405032] . 

A V. Frolov and V. P . Frolov, "Black holes in a compactified spacetime," Phys. Rev. D 67 (2003) 124025 
arXiv:hep-th/0302085 . 

R. Gregory and R. Laflamme, "Black strings and p-branes are unstable," Phys. Rev. Lett. 70 (1993) 2837 
arXiv:hep-th/9301052 . 

R. Gregory and R. Laflamme, "The Instability of charged black strings and p-branes," Nucl. Phys. B 428 (1994) 
399 arXiv:hep-th/9404071 . 

G. T. Horowitz and K . Maeda, "Fate of the black string instability," Phys. Rev. Lett. 87 (2001) 131301 
|arXiv:hep-th/0105111|. 



T. Harmark and N. A. Obers, "Black holes on cylinders," JHEP 0205 (2002) 032 arXiv:hep-th/0204047 . 

T. Harmark and N. A. Obers, "New phase diagram for black holes and strings on cylinders," Class. Quant. Grav. 

21 (2004) 1709 arXiv:hep -th/0309116 . 

T. Harmark and N. A. Obers, "Phase structure of black holes and strings on cylinders," Nucl. Phys. B 684 
(2004) 183 arXiv:hep-th/0309230 . 

T. Harmark "Small black holes on cylinders," Phys. Rev. D 69 (2004) 104015 arXiv:hep-th/0310259 . 

H. Elvang, T. Harmark and N. A. Obers, "Sequences of bubbles and holes: New phases of Kaluza-Klein black 

holes," JHEP 0501 (2005) 003 arXiv:hep-th/0407050 . 

S. S. Gubser, "On non-uniform black branes," Class. Quant. Grav. 19 (2002) 4825 |arXiv:hep-th/0110193|. 
B. Kol, "Topology change in general relativity and the black-hole black-string transition," arXiv:hep-th/0206220 
B. Kol and T. Wiseman, "Evidence that highly non-uniform black strings have a conical waist," Class. Quant. 
Grav. 20 (2003) 3493 arXiv:hep-th/0304070 . 

B. Kol, E. Sorkin and T. Piran, "Caged black holes: Black holes in compactified spacetimes. I: Theory," Phys. 
Rev. D 69 (2004) 064031 arXiv:hep-th/0309190 . 

E. Sorkin, B. Kol and T. Piran, "Caged black holes: Black holes in compactified spacetimes. II: 5d numerical 
implementation," Phys. Rev. D 69 (2004) 064032 arXiv:hep-th/0310096 . 

D. Gorbonos and B. Kol, "A dialogue of multipoles: Matched asymptotic expansion for caged black holes," 
JHEP 0406 (2004) 053 arXiv:hep-th/0406002 . 

D. Gorbonos and B. Kol, "Matched asymptotic expansion for caged black holes: Regularization of the post- 
Newtonian order," arXiv:hep-th/0505009 

T. Wiseman, "Static axisymmetric vacuum solutions and non-uniform black strings," Class. Quant. Grav. 20 
(2003) 1137 |arXiv:hep-th/0209051 



T. Wiseman, "From black strings to black holes," Class. Quant. Grav. 20 (2003) 1177 arXiv:hep-th/0211028 . 
H. Kudoh and T. W iseman, "Properties of Kaluza-Klein black holes," Prog. Theor. Phys. Ill (2004) 475 
|arXiv:hep-th/031uT04] . 

H. Kudoh and T. Wiseman, "Connecting black holes and black strings," Phys. Rev. Lett. 94 (2005) 161102 
arXiv:hep-th/0409111 . 

R. Emparan and H . S. Reall, "Generalized Weyl solutions," Phys. Rev. D 65 (2002) 084025 
arXiv:hep-th/0110258 . 

S. Fairhurst and B. Krishnan, "Distorted black holes with charge," Int. J. Mod. Phys. D 10 (2001) 691 
|arXiv:gr- qc/0010088 . 

R. Geroch and J. B. Hartle, "Distorted black holes," Jour. Math. Phys. 23 (1981) 680 

J. Bicak, T. Ledvinka, B. G. Schmidt and M. Zofka, "Static fluid cylinders and their fields: global solutions," 
Class. Quant. Grav. 21 (2004) 1583 arXiv:gr-qc/0403012 . 
[30] L. Richterek, J. Novotny and J. Horsky, "New Einstein-Maxwell fields of Levi-Civita's type," Czech. J. Phys. 



20 



50 (2000) 925 arXiv:gr-qc/0003004 . 

[31] A. Y. Miguelote, M. F. A. da Silva, A. Wang and N. O. Santos, "Levi-Civita solutions coupled with electromag- 
netic fields," Class. Quant. Grav. 18 (2001) 4569 [arXiv:gr-qc/0104018| . 

[32] M. Nouri-Zonoz, "Cylindrical analogue of NUT space: Spacetime of a line gravomagnetic monopole," Class. 
Quant. Grav. 14 (1997) 3123 arXiv:gr-qc/9706015 . 

[33] S. S. Yazadjiev, "Distorted charged dilaton black holes," Class. Quant. Grav. 18 (2001) 2105 
arXiv:gr-qc/0012009 . 

[34] We use units such that the four-dimensional Einstein-Maxwell equations read G a b = 2T a b where T a b = F a c Fb c — 
\F cd F cc ig a b- 



